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ABSTRACT 

The WR binary CV Serpentis (= WR113, WC8d + 08-9IV) has been a source of 
mystery since it was shown that its atmospheric echpses change with time over decades, 
in addition to its sporadic dust production. The first high-precision time-dependent 
photometric observations obtained with the MOST space telescope in 2009 show two 
consecutive eclipses over the 29d orbit, with varying depths. A subsequent MOST run 
in 2010 showed a seemingly asymmetric eclipse profile. In order to help make sense of 
these observations, parallel optical spectroscopy was obtained from the Mont Megantic 
Observatory (2009, 2010) and from the Dominion Astrophysical Observatory (2009). 
Assuming these depth variations are entirely due to electron scattering in a /3-law wind, 
an unprecedented 62% increase in M is observed over one orbital period. Alternatively, 
no change in mass-loss rate would be required if a relatively small fraction of the 
carbon ions in the wind globally recombined and coaggulated to form carbon dust 
grains. However, it remains a mystery as to how this could occur. There also seems to 
be evidence for the presence of corotating interaction regions (CIR) in the WR wind: 
a CIR-like signature is found in the light curves, implying a potential rotation period 
for the WR star of 1.6 d. Finally, a new circular orbit is derived, along with constraints 
for the wind collision. 

Key words: binaries: eclipsing - stars: mass-loss - stars: winds, outflows - stars: 
Wolf-Rayet. 



1 INTRODUCTION 

Ever since lAllen et all |l9z2) 

showed an IR excess in some 
WC9 stars, it has been known that certain late-type WCs 
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produc e dust. This dust is c omposed of amorphous carbon 
grams l|Williams et all Il987l ') and its formation is also fa- 
vored in suitable WC -I- O binaries (i.e. binaries with large 
enough separations with respect to each star's luminosity, so 
that the ionizing flux from the stars cannot prevent dust for- 
mation), presumably because of the high densities attained 
in the shocked region between the colliding winds. 
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CV Ser (= HD 168206 = WR113, a (J2000.0) = 
18:19:07.36, 5 (J2000.0) = -11:37:59.2, v = 9. 2) is a 
long-studied WC8d-|-08-9IV spectroscopic binary l|Hiltneij 
Il945h with atmosph eric eclipses and a 29.704d period 



i Niemela et al. 19961). Follow ing the first published light- 
curve by iGaposchkinI (| 19491 ) . various other light curves 



whatsoever (e.g. 


Hiellmins: & Hihnei 1961 Stepienl 1970l: 


Kuhi & Schweizer 


Il970l: WiUiams et all Il977l: iLipunoval 


1985h. DifTerent explanations were given, including the dos- 



sibilit y that the authors had used the wrong orbital pe- 
riod jCowlev et al.lll97lh . However, even when no eclipses 
were found in the optical continuum, ICherepashciiu^ (Il972l ) 
showed that the syste m was still e clipsing in the A4653 emis- 
sion line (confirmed by I Morrison fc Wolfi 1972). Since then, 
two MOST (Microvariability and Oscillations of STars; see 
below) runs conducted in 2009 and 2010 also show varying 

depths, possibly implying a varying mass-loss rate. 

CV Ser was also shown to produce dust (jCohen et al.l 

Il97d ). This phenomenon was seen as a plausible explanation 
for the variation in its optical eclipses . It has since been 
classified as a persistent dust producer l|Williamslll995l ). 

Several studies have been carried out to refine the or- 
bital solution of CV Ser. Mo st of them found a quasi circu- 
lar orbit (e.g. iMassev fc Niem cla 1981,) but more recently, 
iNiemela et al.l l| 19961 ) have cast some doubt on that result, 
finding an eccentricity of 0.19. 

Initially, the goal of this project was to monitor stochas- 
tic short-term absorption features in CV Ser's light curve 
involving light from the orbiting O star in order to try and 
link them to clumps in the WR wind. Because of its high 
sensitivity, the MOST space telescope seemed like the per- 
fect instrument to carry out this research. The system was 
observed in 2009 for more than a complete orbital cycle, 
since the anticipated absorption associated with the clumps 
should vary with the orbital phase, depending on the part 
of the WR wind illuminated by the O star along the line of 
sight (see Fig. 1). 

If the clumps produce observable absorption throughout 
the orbit, one would expect to find shorter and deeper vari- 
ations near (f> = Q (WR at inferior conjunction) and almost 
no random variations (unless one of the stars undergoes in- 
trinsic variability) around 4> = 0.5. Therefore, the scatter of 
the residuals from a given light curve fit should vary with 
phase, being larger around = 0. These variations would 
be a few mmag deep, as expected from the typical relative 
density enhancement in clumps. 

Section [2] will briefiy summarize the observations (both 
photometric and spectroscopic) of CV Ser taken for this 
study. The specific results for each observing run are pre- 
sented in section |3] and then are briefiy discussed in sec- 
tion H 



2 OBSERVATIONS 
2.1 Optical photometry 

Two photometry runs (2009 and 2010) were c arried out us- 
ing the MOST space telescope (| Walker et all 12003). a 150- 
mm aperture Rumak-Maksutov telescope with, initially, 2 
CCDs (1 for science and 1 for tracking; however, the track- 
ing CCD was lost in 2007 after being hit by a cosmic ray 



Wavelength (angstroms) 

Figure 2. Sample spectrum from the 2009 OMM dataset. 



and tracking and science images are now both obtained on 
the same CCD) on a polar. Sun- synchronous, low-Earth 101- 
minute orbit. MOST is used for optical photometry, with a 
wideband filter (3000 A, centered around 5250 A). 

The first run was obtained contiguously from 2009 June 
27 to 2009 August 11. This 45-day observation encompassed 
2 eclipses of CV Ser, with millimagnitude precision. Most 
data points are 0.5 min apart, except when the satellite went 
over the South Atlantic Anomaly (SAA). Part of the run 
was also shared with another target, giving an overall some- 
what inhomogeneously distributed data set. One particular- 
ity of this light curve is the difference of depth between both 
eclipses, which will be discussed in more detail in the next 
section. 

Because of this interesting behavior, CV Ser was ob- 
served once again with MOST the following year, from 2010 
June 14 to 2010 July 11, in similar conditions. However, this 
run covered only one eclipse of CV Ser. 



2.2 Optical spectroscopy 

Spectroscopy was obtained at both Observatoire du Mont- 
Megantic (2009 and 2010) and the Dominion Astrophysical 
Observatory (2009 only). Details of these observations are 
summarized in Table [1] 

Unfortunately, some of the 2010 OMM spectra are full 
of narrow instrumental defects. As this effect could not be 
corrected, these spectra had to be discarded. The spectra 
were reduced using IRAF 0. All spectra were normalized to 
the continuum. An example spectrum from OMM is shown 
in Fig. 2, revealing the main WC8 spectral features. 



3 ANALYSIS AND RESULTS 
3.1 MOST 2009 Dataset 

First, without loss of information, the data were binned for 
each MOST orbit (101 minutes), applying 3-sigma clipping 
for each binned point. The resulting light curve is shown in 
Fig. 3a. A steady linear drift can be seen throughout the 
light curve (almost 0.01 mag over 45 days). It did not seem 



^ IRAF is distributed by tiie National Optical Astronomy Obser- 
vatories (NOAO), which is operated by the Association of Univer- 
sities for Research in Astronomy, Inc., (AURA) under cooperative 
agreement with the National Science Foundation (NSF). 
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Phase =0 




Phase = 0.5 



Figure 1. Intuitive cartoon to show how it is expected that the absorption by the clumps along the line of sight should depend on the 
orbital phase. 



Table 1. Summary of the spectroscopic observations. 



Telescope 


Year 


Dates 




Number of spectra 


Spectral range (A) 


Resolution (A/pixel) 


Average SNR 


OMM (1.6m) 


2009 


July 5 - 


August 8 


145 


4400-6000 


~ 0.64 


~ 150 


DAO (1.8m) 


2009 


July 7 - 


July 31 


156 


5300-6100 


~ 0.77 


> 100 


OMM (1.6m) 


2010 


June 21 


- July 12 


50 


4400-6000 


~ 0.64 


~ 150 



to be due only to the stars and it is well known that such 
instrumental drifts occur over long periods of continuous ob- 
servation with the MOST satellite. Similar drifts are visible 
in the light curves of other field stars, but since the slope 
of these drifts does not vary smoothly from one region to 
the next on the CCD, it is impossible to map them out in 
an effort to determine what part of the drift of the CV Ser 
light curve is really instrumental. The only option was to fit 
it out. Nevertheless, both eclipses had undeniably different 
depths, unless an unusually deep stochastic dip (e.g. due to 
a large clump), like those seen elsewhere in the light curve, 
happened to occur exactly at phase 0.00. Not only does such 
a coincidence seem extremely unlikely, but a stochastic dip 
would not likely produce the same eclipse profile due to the 
global WR wind. Therefore, we dismiss this possibility as 
being implausible. 

The data were fitt ed using the model developed in 
iLamontagne et akl (|l996l ) : the eclipse is atmospheric and oc- 
curs as the O star goes behind the WR wind, its light being 
Thomson scattered by free electrons in the WR wind. This 
assumption is reasonable, because at the considered wave- 
lengths, Thomson scattering is clearly the dominating pro- 
cess. Compared, for example, to free-free absorption, Thom- 
son scattering's infiuence on the light curve is more than one 
order of magnitude greater in the visible spectrum, even in 
the densest wind regions (at R = 2R,). Only in the NIR do 
these processes become comparable. This is easily verifiable 



using the appropriate formulae of iTuckerl |l973)- A circular 
orbit is also assumed in this model. 

The Lamontagne et al. model is purely geometrical and 
it yields various parameters of the system, including the or- 
bital inclination and the mass-loss rate of the WR compo- 
nent. One of the advantages of this model is that the ob- 
tained mass-loss rate is not affected by clumping, and since 
it makes simple assumptions, it is fairly easy to use. 

For wind velocity laws with integer values of /3, the 
model produces an analytical formula describing the behav- 
ior of the light curve. In this case we used a /3 = law (i.e. 
constant expansion velocity), although it is not the closest 
to reality (at least not in the case of hot stars, where we 
have P ^ 0.8, Pauldrach et al .111 9861 ). However, as shown in 



ILamontagne et al] (| 19961 ). both P = and 13 — 1 velocity 
laws in practice produce very similar synthetic light curves, 
and the first case requires fewer free parameters. The corre- 
sponding formula is 

A A ^ f 7r/2 + arcsine\ 

Am = Amo + A ' — (1) 



with Amo a constant, e = (sini) cos27r(;/), 
A = J^^iHi^ and fc - "'^^^^ 



( 1 -I- IwR / lo ) Anm-pVoo a ' 

in which i is the orbital inclination, e is Euler's number 
(~ 2.718), Iwr/Io is the intensity ratio of the 2 stars in the 
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Table 2. Best fit values for the modified Lamontagne model based 
on the 2009 CV Ser light curve. 



Model parameter Best fit value 



Bo (maf 


;) 


-0.0101 ± 0.0004 


Bi (maf 


;/P) 


0.0040 ± 0.0005 


^0 (maf 


;) 


0.0013 ±0.0001 


Al (maf 


;/P) 


0.0008 ± 0.0001 


sini 




0.979 ± 0.003 



observed bandpass, a ~ 0.5 free electron per baryon mass 
(since the dominating element is fully ionized helium), cje is 
the Thomson electron-scattering cross section, M is the WR 
mass-loss rate, m-p is the proton mass, Voa is the terminal 
wind speed and a is the orbital separation. 

However, there was one small problem trying to apply 
this model to our data, since it only accounts for constant 
winds, whereas the changing depth of the eclipses in CV 
Ser's light curve clearly indicates that somehow, over one 
orbital cycle, the parameters of the WR wind have changed. 
In order to take this effect into account, it was necessary 
to slightly modify the model to include possible variations 
to the A parameter, which includes most of the important 
parameters of the WR wind. Since little is known about 
the processes involved in such a change in the wind's struc- 
ture, there was no reason to assume a particular shape for 
the variation of A, so we only included a first-order, time- 
dependent term (which also seems justified by the relatively 
short period of time covered by the data) . 

Because of the instrumental drift, we also allowed Amo 
(replaced with Bq in the formula) to vary linearly with time 
(with slope Bi), giving the following equation for the mod- 
ified model: 

„ „ / . .si 7r/2 -f- arcsin e\ , , 

/\m = Bo + Bit + [Ao + Alt) \-^-j== — j, (2) 

where t is given in units of phase ([1/P]). 

The corresponding fit is plotted over the light curve in 
Fig. 3a. The values of the different parameters are given 
in Table [2l A Monte-Carlo simulation was used, randomly 
distributing the errors on the data, in order to evaluate the 
uncertainties on the fit parameters. 

One thing particularly stands out when looking at these 
values: the ratio A\P/Ao ~ 0.62 ± 0.12, which means that 
the A parameter, which was considered to be constant in 
the initial model, incre ases by about 6 2% over one complete 
orbital cycle (29.704d. iNiemela et"aLlil996 '). But how does 
this afTect the values of the physical parameters of the WR 
wind? We must first look at each parameter contained in A 
individually: the ionization level (a), the terminal velocity 
(woo), the intensity ratio (Iwr/Io) and the mass-loss rate 
(M). 

The first three parameters all should affect the spectra 
of the system. However, as it will be discussed in subsec- 
tion [3]3] no such changes were detected. 

Therefore, we conclude that only the mass-loss rate 
varies significantly and its variation is of the order of 62% 
over about 30 days, which is unprecedented in any WR 
star. Using the same (constant) values of a = 0.5, Voo = 
1890k m/s and Iwr/Io ~ 0.69 as in iLamontagne et al] 
l|l996t ). we obtain a value of M = 3.5 x IO'^Mqi/"^ for 



the first eclipse, and Ad = 5.7 x 10~^ MQy~^ for the second 
eclipse. 

Next, the fit was subtracted from the light curve and 
the residuals were analysed to see if there were any signs 
of random variations due to clumping. No significant phase- 
dependent change in the level of random variability was de- 
tected, making a link to clumping unlikely. 

The Fourier an alysis of the residuals (using period04, 
iLenz fc Bregeill2005l ) yielded a significant peak with a 1.18± 
0.03 c/d frequency (see Fig. 4a). This frequency was ruled 
out as an instrumental effect and therefore seemed to be 
intrinsic to the system. Time- frequency analysis (Fig. 5a) 
provided a clearer picture: both a ~ 1.18 c/d frequency 
and a broad ~ 0.6 c/d frequency appeared intermittently 
(and alternately). This signature was very si milar to that 
detected in the time- frequency plot of WRllO (jChene et al.l 
l201lh . suggesting the presence of corotating interaction re- 
gions (CIRs) in WR113's wind. 

The obtained frequencies seem to be coherent with the 
hypothesis th at CIRs may com e in pairs (or even higher 
multiples, e.g. iKaper et al.lll997l ') per rotation cycle. Indeed, 
since one of the frequencies is about twice the other, as is the 
case for WRllO, it seems plausible that the lower frequency 
is associated with the star's rotation period, whereas the 
higher one is due to the presence of two CIR arms. However, 
in order to be able to conclude that the star rotates at a 
frequency of ~ 0.6 c/d, it would be preferable to see if this 
frequency is steady over longer periods of time. 

3.2 MOST 2010 Dataset 

The 2010 MOST light curve was binned exactly like that 
from 2009, in order to be able to make an easier comparison 
between both light curves. It consists of 27 days of continu- 
ous coverage, spanning across one eclipse (Fig. 3b). 

The eclipse shown in this light curve seems to be 
asymmetrical, a shape some what reminiscent of the light 
curve published by HicUmi ng fc Hiltnej (119631 ). although it 
is much shallower. The same extra absorption during ingress 
is present. One hypothesis could be that this additional ab- 
sorption might be due to a dust event. However, it isn't clear 
whether this depression is actually statistically significant or 
not, since it seems to be quite weak. In order to test its sig- 
nificance, the standard Lamontagne et al. model was used 
to fit the eclipse, once using the entire light curve, and once 
using only the second half of the eclipse to avoid including 
this depression. Both fits were done using MPFIT (a robust 
non-linear least squares curve fit ting routine based on the 
Levenberg-Marquardt algorithm, iMarkwardtl 120091 ) and in 
both cases, the value of sin i was fixed to that found in 
2009. 

The original model (i.e. with constant A parameter and 
linear instrumental drift) was used instead of the modified 
one presented earlier in this article since there was only one 
eclipse in the 2010 light curve. Indeed, outside of the eclipses, 
a linear variation of the A parameter induced an effect simi- 
lar to a simple linear drift, since both parameters are corre- 
lated (this was made obvious by the Monte-Carlo simulation 
used to determine the errors on the fit parameters for the 
previous dataset), so in the case of the 2009 light curve, the 
only thing making it possible to constrain the value of that 
variation was the presence of 2 eclipses. Therefore, this im- 
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Modified Lamontagne fit - WR113 



Lamontagne fit - 2010 MOST iight curve 




0,04 i 

-0.4 



(a) Modified Lamontagne fit to the 2009 MOST light curve. 



(b) Lamontagne fit to the 2010 data. 



Figure 3. MOST hght curve of CV Ser (a ) 2009: the most str i king f eature is the difference between the depth of each eclipse. Overplotted 
is the fit to the data, using the modified iLamontagne et al. model, as explained in the text. A typical single-point 2cr errorbar 

(also valid for 2010) is shown, (b) 2010: the eclipse seems to be asymmetrical. Overplotted is the fit to the data (as described in the 
text), using the original Lamontagne model with a linear drift term. The phases were computed using a 29.700d period. 




3.0 4.0 
Frequency 



(a) 2009 periodogram. 



0.0012 



j 0.0008 



0.0004 



O.OOOi 




(b) 2010 periodogram. 



Figure 4. Periodogram of the 2009 MOST light curve residuals (a): a peak can be found at a frequency of 1.18 c/d. Periodogram of the 
2010 MOST light curve residuals (b): the standout peak is at a frequency of 0.64 c/d. (The amplitude is in mag and the frequency in 
c/d for both plots.) 



possibility of disentangling the linear drift and the variation 
of A renders any attempt of evaluating Ai futile. A can only 
be found accurately during the eclipse, and the value of the 
instrumental linear drift is not relevant (since it cannot be 
realistically constrained). 

Within the errors, both fits yielded the same values for 
the model parameters, thus suggesting that the asymmetry 
is dominated by a random artefact rather than a truly sig- 
nificant feature of the hght curve. The best fit values of the 
standard Lamontagne model (including a parameter to de- 
scribe once again the instrumental linear drift) are shown in 
Tabled 

Deriving the mass-loss rate during the eclipse from the 
fitted value of the A parameter, as for the 2009 data, one 
gets a value of M = 2.8 x 1Q~^ Mqu'^ , which is the value of 
the mass-loss rate during the eclipse. This value is slightly 



Table 3. Best fit values for the Lamontagne model based on the 
2010 CV Ser light curve. 



Model parameter 


Best fit value 


Bo (mag) 


-0.0040 ± 0.0003 


Bi (mag/P) 


0.0058 ± 0.0006 


A (constant, mag) 


0.00104 ± 0.00005 


sini (fixed) 


0.979 



inferior to the mass-loss rate observed during the first eclipse 
of 2009 (and quite inferior to the rate observed during the 
second eclipse), suggesting that this quantity might vary 
considera bly with time in the case of CV Ser (a s seen in 
iHiellming fc Hiltnei|[T963 : iKuhi fc Schweizeilll970l . etc.). 
The fit was then subtracted from the hght curve in an 
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(a) 2009 time-frequency plot. (b) 2010 time-frequency plot. 

Figure 5. Time-frequency plot of the MOST light curve residuals (a)2009: a CIR-like signature is detected. (b)2010: a similar signature 
to that in 2009 is found, strengthening the CIR hypothesis. 



Lomontaqne fit residuals - 2010 MOST light curve 



Figure 6. 2010 MOST light curve residuals. The arrows corre- 
spond approximately to the inferred rotation period. Black arrows 
show peaks in the residuals, possibly linked to CIRs, whereas no 
such peaks are found where there are grey arrows (it should be 
recalled that the CIR phenomenon is expected to be cyclical and 
not strictly periodic). Similar peaks are found in the 2009 resid- 
uals. 



effort to try to isolate the cause of the asymmetry (possibly a 
dust event), but it is not clear from the residuals plot (Fig. 6) 
that the depression is actually anything but noise (because 
it is much shallower than the surrounding variability) . 

As in the case of the 2009 hght curve, no significant 
signs of random phase-dependent variations were found. 
The Fourier analysis of the residuals (Fig. 4b) yielded a 



0.64 ± 0.03 c/d peak, coherent with one of the frequencies 
found in 2009. The time-frequency plot (Fig. 5b) shows the 
same type of signature as the 2009 hght curve: two main 
frequencies (around 0.64 c/d and 1.1 c/d) come and go dur- 
ing the 27-day time series. The same frequencies were found 
the previous year, which is coherent with the CIR scenario, 
since the apparent frequencies should be modulated by the 
rotation rate. Therefore, this suggests a rotation period of 
1.6d for the WR star according to the argument previously 
stated, thus giving a r atio of Vrot/vcrit ~ 22%, which is rela- 
tively high for a WR, iHarries et al.l 19981 (vcrit is calculated 
using 7?, — 5Rq, as in iLamontagne et al.l[l99^ 1. It should 
also be mentioned that a lower frequency is found approx- 
imately around the time of the apparent depression in the 
light curve (lower left of Fig. 5b). It lasts for about 6 days 
and corresponds to a period of about 3-4 days. However, it 
is hard to determine whether it is linked to the depression 
or not. 



3.3 Spectroscopy 

One of the goals of this spectroscopic campaign was to 
solve the orbit and det ermine wheth er it is circular (e.g . 
iMassev fc Niemela|[l98ll ) or eccentric l|Niemela et al. mm . 
The radial velocities were obtained using the FXCOR rou- 
tine in IRAF over the entire range of the DAO spectra, since 
tests using different ranges show that the excess emission in 
Cm A5696 and the O star's absorption lines do not perturb 
significantly the obtained WR RVs. The same range was 
used with the OMM spectra for consistency. Fig. 7 shows 
the relative RV plot, as well as the best circular fit. Unfor- 
tunately, the O star's absorption features being negligible 
and very weak in comparison to the strong, broad emission 
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Figure 7. Relative radial velocities of the WR component ac- 
cording to phase. The fitted circular orbit is overplotted. The 
dots represent the DAO spectra while the asterisks represent the 
OMM data (black for 2009, grey for 2010). 



Figure 8. Equivalent width of the Cm A4650 (top) and Hell 
A4686 (bottom) lines as a function of phase, based on OMM (2009 
= black, 2010 = grey) spectra. A double gaussian fit was applied 
in order to deblend both lines, resulting in some minor artefacts. 



Table 4. Best fit values for the WR orbit based on the current 
data only. 



Orbital parameter 


Elliptical fit 


Circular fit 


To (HJD) 


2455003 ± 1 




E (HJD) 




2455012.6 ±0.1 


K (km/s) 


173 ± 1 


172 ± 1 


e 


0.02 ±0.02 


(fixed) 


^(°) 


330 ± 10 




P(d) 


29.704 (fixed) 


29.700 ±0.001 



lines of the WR star, it was not possible to establish that 
star's orbit and Fig. 7 only shows the WR orbit. 

A first elliptica l fit w as done using a 29.704d period 
as in iNiemela et al.] (|l996l ). However, unlike in that article, 
we obtain an orbit which is indistinguishable from circu- 
lar. Hence, a circular fit was then applied, and in conjunc- 
tion with data from iMassev fc Niemelal (|l98ll ). we obtain 
a revised period of 29.700d. The obtained value of K is 
also somewhat larger than previously report ed values (e.g. 
iMassev fc Niemela|[l98ll : iNiemela et al.ll 19961 ) . The parame- 
ters of both fits are shown in Table [l] 

A closer look at the individual spectra reveals consider- 
able spectral variability. There are different sources of vari- 
ability which can be distinguished from one another. First, 
there is short-term variability due to the clumped material 
in the WR wind (including the wind collision zone) . This is 
seen as small, variable subpeaks atop some emission lines. 
Then, some emission lines (Cm A5696 in particular) show 
excess emission coming from the wind collision zone between 
the WR and O components. As the shock cone orbits with 
the stars, the excess shifts from one side to the other of the 
emission line, with a quarter-phase delay. This effect will be 
studied further below. 

Finally, it is well kn own that some emission lines eclips e 
while others do not (e.g. lCherepashchukll97j : ICowlevll 19721 ). 
This behavior is exhibited in the spectra obtained for this 
study, for instance in the case of Cm A4650 (as illustrated in 



Fig. 8, this line eclipses while another line, Hell A4686, does 
not). Even though some lines eclipse, it appears obvious that 
there is no dramatic increase or decrease in their EW outside 
of these eclipses (as shown in Fig. 8 for Cm A4650). There- 
fore, the ratio of the intensities of the O star and the WR 
star must be quite constant, since e.g. if Iwr/Io were to 
increase, the emission lines of the WR star would be diluted 
by the O star's continuum and would have lower equiva- 
lent widths. Thus, we conclude that Iwr/Io does not con- 
tribute significantly to the variation of the A parameter, as 
discussed in section [3711 

Another wind parameter that was monitored is the ter- 
minal velocity, fn order to trace any possible variation of 
«oo, it is best to use a low-ionization line formed far out in 
the wind. In this case, Hel A5876 was used, but no signifi- 
cant variation was found. Thus, Voo cannot account for the 
variation of the A parameter either. 

Finally, following the study of the equivalent widths of 
individual lines, one can check the ratio between lines of dif- 
ferent ionization levels, to see if the spectral type changes 
with time. Two important lines in the determination of 
WC subtypes are Cm A5696 and Civ A5808. The ratio of 
their equivalent widths (including the excess emission in Cm 
A5696 coming from the wind collision zone, which has con- 
stant flux) is fairly constant, implying a constant spectral 
type and, therefore, ionization level. Consequently, as dis- 
cussed previously, we conclude that the only wind param- 
eter that plays a role in the variation of the depth of the 
eclipses in 2009 is the mass-loss rate. 

As mentioned previously, excess emission from the 
wind-collision zone may o ffer a lot of information. In this 
case, we used the model of iLiihrsI l|l997h to study the excess 
emission in the Cm 5696 line. 

The Liihrs model considers the wind collision zone, 
which takes the form of a shock cone (with a rounded head) . 
This corresponds to the surface where the winds of both 
stars collide with equal and opposite momentum compo- 
nents (although more emission may come from the WR side 
of the shock) . According to the behavior of the excess emis- 
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sion, it is possible to determine 9, the half-opening angle of 
the cone, 59, the thickness of the cone, 54>, the phase delay 
induced by the Coriolis force acting upon the shock cone as 
both stars orbit, Ustrm, the streaming velocity of the mate- 
rial as it flows along the cone and i, the orbital inclination. 

While the Liihrs model requires an intricate fit of each 
individu al spectrum, ther e also exists an integral form, as 
found in iHill et all (|200(]| '). This form only requires the RV 
and FWHM (full width at half maximum) of the excess emis- 
sion as a function of phase. The two main equations are the 
following : 

FWHMe.:c = Ci+2vstrm siu g ^1 - siu^ j cos2 (0 - 50) (3) 
and 

RVex ~ C2 + Vstrm COS 9 siu j COS {(f> — S(f>) . (4) 

This form of the model does not yield the thickness of 
the cone (59). 

In addition to studying the wind collision zone, which 
in itself is an interesting phenomenon, one might be tempted 
to try to find proof for the variation of Mwr by applying 
this model, since there exists a relationship between the half- 
opening a ngle o f the s hock cone and the ratio of the wind 
momenta. lUsovl ([l99^ finds the following formula : 

9irad)c^ 2.1 (^1 - liPj r,'^\ (5) 

where n — , ^o^x,,o — comprised between 10~* and 

1. If the half-opening angle of the shock cone can accurately 
map out the variations of the WR's mass-loss rate, such a 
spectral analysis could prove useful in order to monitor a 
system like CV Ser, where M seems to vary quite a bit. 

Fig. 9 shows the Cm A5696 line in the WR reference 
frame. It is quite obvious here that the overlying excess 
varies, but in order to be able to isolate it, a minimum profile 
must be subtracted, which should correspond to the unper- 
turbed emission line. Here, the line seems a little slanted, so 
the base profile (shown in Fig. 9) was chosen to be a tilted 
flat top (in order to fit both the wings and the unperturbed 
edges of the top of the line). Once subtracted, it is easier 
to see how the excess varies with phase (Fig. 10). 

The Liihrs analysis was first performed by calculating 
the first moment (mean) of the excess for the RV and roughly 
twice the square root of the second centered moment (stan- 
dard deviation) for the FWHM. One difficulty encountered 
using this technique is that the result varied greatly depend- 
ing upon what integration bounds were used. Consequently, 
this technique was not reliable enough, using the available 
data, and had to be abandoned and replaced by a simpler 
one - actually measuring the FWHM and using the midpoint 
at half maximum to compute the radial velocities. 

To do so, the OMM and DAO spectra were combined 
into a single dataset, which was then binned to get rid of 
short-term stochastic effects as well as diminish the noise. 
Each individual bin covers 1/20 of the orbital period. The 
inclination angle, i, was fixed to the value found using the 

^ Although the choice of such a profile may seem a bit ad hoc, 
similar profiles hav e been found in real sp ectra (e.g. in the single 
WC8 star WR135. iLipine fc MoffagilQQgh . 
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5660 5680 5700 5720 5740 
Wovelength (ongstroms) 

Figure 10. Excess emission in Cm A5696 varying with phase 
along the vertical axis. Combined data from OMM (2009, 2010) 
and DAO (2009), binned by phase (bins of 0.05). 



photometry, since it is not well constrained by Eqs. [3]and|3] 
in the presence of instrumental noise. The results are shown 
in Fig. 11. 

The fit was first performed without letting 9 vary with 
time. The results are presented in Table [S] Unfortunately, 
since the dependence of S on 77 is rather weak, the slight 
effect it would have on the data is too subtle to be well 
constrained; therefore it was not possible to highlight the 
effects of the varying mass-loss rate on the shock cone. Also, 
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CHI 5696 base profile 




L , , , \ , , , \ , , , \ , , , A 

5660 5680 5700 5720 5740 

Wavelength (angstroms) 



Figure 9. Base profile (grey line), which is to be subtracted from the Cm A5696 line to isolate the excess emission. The Cm A5696 
profiles from both the DAO (2009) and OMM (2009, 2010) data are shown, shifted in the WR reference frame, with their base profile. 
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Figure 11. FW (top) and RV (bottom) of the excess emission 
as a function of phase. Each dot correspond to a particular phase 
bin (the OMM and DAO data were combined together for this 
analysis) . 



the obtained values are somewhat surprisi nRly low and far 
from those found in the prev ious study of 'Antokhin et al. 

Using the lUsovl (|l995l ) formu la, an rj param eter of 
about 0.012 is found, whereas using iMartins et all (|2005l ) 
(to obtain reasonable stel lar parameters acc ording to the O 
star's spectral type) and IVink et al.l l|200ll ) (to obtain Mo 
and v^^o in function of those stellar parameters), we expect 
to find a value of about 0.08. The two obtained values are 
a little less than an order of magnitude apart, which seems 



Table 5. Best fit values for the Liihrs analysis (based on the DAO 
spectra). 



Model parameter 


Best fit value 


Vstrm (km/s) 


366 ± 18 


6 (deg) 


27 lb 3 


S<f> (deg) 


36 ±2 



fairly reasonable given the many possible sources of error for 
such a calculation. 

Finally, one of the secondary aims of the spectroscopic 
part of this study was to confirm the spectral types of 
both orbital companions. In order to disentangle the spec- 
tra of both stars, th e "shift and add" method described in 
iDemers et al] (|2002l ) was to be used. This method is iterative 
and uses the orbits to "shift" the spectra to the reference 
frame of the WR star. Then, the spectra are "co-added" and 
the resulting compound spectrum is then shifted back to our 
reference frame to be subtracted from the initial spectra, re- 
vealing the spectra of the O star. After a few iterations, it 
is possible to get the spectrum of each star. Since the or- 
bit of the O star was not obtained in this study, we used 
the orbit found in a prev ious study with similar ephemeris 
dMassev fc Niemela|[l98ll ). The method was applied over 10 
iterations. However, the obtained O spectrum is heavily con- 
taminated by artefacts due to the strong WR emission lines. 
This may be partly due to a slightly wrong orbital solution. 
In any case, it is impossible to confirm the spectral classifi- 
cation of the O companion based on these results. 

However, the WR spectrum does confirm the previous 
WC8d classification (the 'd' denotes the detection of dust 
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in the IR), since it corresponds to the features found in 
ISmith et all (|l990l ). In particular, the Civ A5808/Ciii A5696 
peak flux ratio is about 0.7, which lies well within the 0.5- 
2.0 range associated with the WC8 spectral type (this ratio 
goes up to ~ 0.8 if the excess emission produced in the wind 
collision zone is removed). 

It should also be mentioned that no attempt was made 
to detect the hypothetical CIRs in the spectroscopic data, as 
the signal is complicated by the presence of excess emission 
from the wind collision zone. 



4 DISCUSSION AND CONCLUSIONS 

After having eliminated the other possible causes for the 
change in eclipse depth observed in the 2009 MOST photom- 
etry, the deduced 62% increase in M over one 29.704d or- 
bital period is a truly remarkable finding, unprecedented for 
WR stars. Even more intriguing is the fact that the follow- 
ing year, the mass-loss rate had gone back to a lower value, 
suggesting that it may vary considerably over short to long 
timescales. The derived values of M are also f airly low for a 
WR star. According to Crowther et all l|2002l '). variations in 
wind density should cause the spectral type to change, but 
this effect is not observed here, adding to the mystery, al- 
though these results are not necessarily incompatible since it 
is noted that the dependence on the wind density is weaker 
for late- type WC stars (WC8 and WC9). In any case, this 
definitely constitutes a challenge to both theorists and ob- 
servers, since there is a dire need for a theoretical expla- 
nation as to what could drive such an impressive variation 
of the mass-loss rate, while constant monitoring would be 
required in order to get a better idea of the long-term be- 
haviour of this parameter. One alternate explanation could 
conceivably be the presence of thin-shell instabilities in the 
wind collision zone. However, this will only affect the light 
curve dips stochastically, not systematically. In any case, it 
is not so much the details of the wind collision which deter- 
mine the eclipse's shape and depth; rather it is the global 
wind of the WR component in which the O star orbits. 

As for the primary goal of this study, it was not possible 
to systematically link the stochastic photometric variability 
to the clumping phenomenon. There was no clear phase de- 
pendence of the scatter of the light curves. The variability 
is then probably due to noise (the combination of the un- 
quantifiable absorption due to clumping and other stochas- 
tic phenomena), but also to CIRs in the WR wind, as the 
analysis of the light curve residuals suggests. The repeated 
signature (both in 2009 and 2010) offers a rather robust clue 
towards that hypothesis. This finding is very interesting in 
the actual context of CIRs in WR stars, since more and more 

St-Louis et all |2009|. 



CIR candidates a. r e showing up (e.g 



IChene et al.l I2OIII . IChene fc St-Louij 



20111 ). Even though 



they have mainly been considered as an exception up un- 
til now, CIRs in WR winds might just prove to be the rule 
after all. Further studies will be necessary in order to find 
cyclical spectroscopic evidence supporting this possibility. 
However, their origin remains mysterious, especially since 
no magnetic fields have been detected in the extended atmo- 
spheres of WR stars yet. Nevertheless, this scenario should 
not be ruled out since at the base of the winds, where the 
CIRs originate, the magnetic field might have much higher 



values than our current detection limits but still is not de- 
tected simply because this region is obscured by the dense 
wind (or perhaps due to cancellation in small dipole loops as 
seen on sunspots). Non radial pulsations (NRP) or starspots 
have also been suggested as possib le causes for this phe- 
nomenon (|Cranmer fc Owockilll99^ ) . 

Finally, what was initially thought to be a possible dust 
event during the 2010 observations cannot conclusively be 
determined as statistically significant. However, the per- 
ceived asymmetry in the light curve is reminiscent of the 
model presented by I Veen et al.l l|l998h . Although this model 
was used to characterize the condensation of dust clouds 
and its effect on a light curve, and as such is very unlikely 
linked to these observations (since such a dust event should 
create a much deeper dip in the light curve and the proba- 
bility of it occurring exactly at the same time as the eclipse 
is rather low), it still gives way to a very important question 
: what role can dust play in CV Ser's bizarre behaviour? 
Could a variation of the quantity of dust persistently pro- 
duced by the WR star affect the light curve in the same way 
it has been inferred that a variation of the total mass-loss 
rate might? In order for such an explanation to arise, one 
must first postulate that the dust distribution in the wind 
is the same as the free electron and hence overall density 
distribution, in order to preserve the system's geometrical 
properties and to reproduce the same eclipse profile found in 
the light curve. As to whether this assumption makes phys- 
ical sense, it is hard to determine whether the dust could 
be produced (or maintained) following such a spatial distri- 
bution, since the production of dust in late WC stars is not 
yet well understood. However, it is fairly straightforward to 
calculate what quantity of dust would be needed to produce 
a similar effect to that of an increased mass-loss rate on the 
light curve. Assuming Mie scattering, we get a cross-section, 
for one grain of average radius a, of ad = Qna^, where Q 
is the cross-se ction efficiency. Choosin g a to be ~ 0.1 fim 
(as in WR140: iMarchenko et aL I l2003t ). we find that g ~ 2 
for wavelengths of the order of 5000 A, giving a value of 
ad ~ 6.3 X lO^^'' cm^. However, using a density inside the 
dust grains oi pg — 2 g/cm"^, we find that each grain is 
composed of roughly 4 x 10* carbon atoms. Since the most 
common carbon ion in the WR wind is Cm, for each atom 
of carbon we should expect to find 2 free electrons. Then, 
the Thomson cross-section of the number of free electrons 
needed to combine with 4 x 10* carbon atoms to produce 
one grain is ae,tot = 8 x 10* x tJe = 5.3 x 10"^'' cm^. There- 
fore, scattering by the equivalent quantity of recombined 
carbon in the form of a dust grain is about lO'' times more 
efficient, therefore an appreciable increase in the depth of 
the eclipse could be caused by recombination of a negligibly 
small quantity of carbon and its condensation in the form 
of dust. Possible support for this is seen in the UBV light 
curve of CV Ser obtained from 1984 to 1994 by Dzhapi- 
ashvili (Anthokhin, priv. com.) in which the eclipse depth 
varies strongly with wavelength. Therefore, if dust plays any 
role in these varying eclipse depths, the 62% value for the 
increase of M over an orbital period in the 2009 data ob- 
tained in this study should probably be considered as an 
upper limit. Unfortunately, not much more can be deduced 
about the production of dust in CV Ser, except that it re- 
mains a very interesting phenomenon and should be further 
studied. Indeed, with its troubled history and intriguing be- 
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haviour, CV Ser might prove to be a key system for under- 
standing the production of dust in WC+O binaries. It does 
not appear clear whether this process is necessarily due to 
the wind collision zone or if it originates in wind shocks and 
is then intrinsic to the WR component. 

In conclusion, this work possibly raises more questions 
than it answers, but we conclude without a doubt that CV 
Ser is a very important system that might hold the answer 
to old problems. Hopefully, our findings will motivate the 
community to take a deeper look into this remarkable object 
in the years to come. 
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